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ABSTRACT 

The simultaneous and coupled evolution of horizontal branch oscillation (HBO) and normal branch 
oscillation (NBO) in Z-type sources suggests that the production of HBO is connected to NBO and 
is caused by changes in the physical/radiative properties of the inner accretion disk, although there 
is a lack of substantial spectral evidence to support this. In this Letter, we present the results of an 
analysis of a RXTE observation of a Z source GX 5-1, where the 6 Hz NBO is simultaneously detected 
along with a HBO at 51 Hz. The variations in the intensity and the associated power density spectrum 
indicate that the HBO and NBO are strongly coupled, originating from the same location in the inner 
accretion disk. The absence of HBO and NBO in the lower energy bands, an increase in the rms 
amplitude with energy and a smooth transition among them suggest that they are produced in the 
hot inner regions of the accretion disk. Based on a spectral analysis, we found a signature of changing 
or physically modified inner disk front during the coupled HBO and NBO evolution. We explore the 
various models to explain the observed phenomenon and propose that the NBO is affiliated to the 
oscillations in the thick/puffed-up inner region of the accretion disk. 

Subject headings: accretion, accretion disks - binaries: close - stars: individual (GX 5-1) - X-rays: 
binaries 



1. INTRODUCTION 

Various mechanisms have been suggested to explain 
the quasi-periodic oscillations (QPOs) in neutron star 
low mass X-ray binaries (NS LMXBs), but the location 
of their production region in the accretion disk and the 
connecting physical/radiative processes are not yet clar- 
ified. Of the two broad classification of NS LXMBs 
viz. atoll and Z-type sources, the latter has a high 
X-ray luminosity of ~ 10'^® erg s~^ (close to the Ed- 
dington limit), which is secularly varying along three 
main branches forming a Z shape in the hardness inten- 
sity and color-color diagram. The associated branches 
from top to bottom are horizontal branch (HB), normal 
branch (NB) and flaring branch (FB) and NB connects 
the other two (Hasinger & van der Klis 1989; Hasinger 
et al. 1990). Broad features in the power density spectra 
known as quasi-periodic oscillations (QPOs) are often ob- 
served along the branches. Apart from twin kHz QPOs 
(van der Klis 1998), different low and intermediate fre- 
quency QPOs have been observed in Z-type NS LMXBs 
named after their branches viz. horizontal branch oscil- 
lations (HBOs), normal branch oscillations (NBOs) and 
flaring branch oscillations (FBOs) (van der Klis 2006). 
The systematic Z-track evolution and the different os- 
cillations are considered to be caused by changes in the 
mass accretion rate but a few studies seem to rule out any 
relation to the accretion rate (Lin et al. 2009; Homan et 
al. 2010). 

As the Z sources traverse from HB to NB, the HBO 
smoothly changes to NBO (typically a ^50 Hz peak is 
converted into a ^6 Hz peak), though the physical causes 
in the accretion disk, responsible for this variation, are 
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unknown. There are two main models which generally 
explain the observed features of NBO. In the radiation- 
hydrodynamical model (Fortner et al. 1989; Miller & 
Lamb 1992), NBO are associated with the oscillations in 
the optical depth of radially inflowing hot material. In 
the other model, NBO is caused as a result of Keplerian 
rotation of thick accretion disk, with variations in its 
dynamical properties (Alpar et al. 1992). 

Many models were proposed to explain the generation 
of HBOs. In the magnetospheric beat-frequency model, 
HBO arises due to the orbital motion at a suitable mag- 
netospheric radius (Alpar & Shamam 1985; Lamb et al. 
1985). The sonic-point beat frequency model, an exten- 
sion of the magnetospheric beat-frequency model, was 
proposed to explain the simultaneous presence of HBOs 
and kHz QPOs (Miller et al. 1996, 1998). In the rela- 
tivistic precession model (Stella & Vietri 1999; Stella et 
al. 1999), HBOs and kHz QPOs are caused by the rel- 
ativistic frame dragging and relativistic periastron pre- 
cession of an eccentric ring. In the two oscillator model 
(Titarchuk & Osherovich 1999; Titarchuk et al. 1999), 
HBO and kHz QPOs are the result of oscillations in a 
hot blob and the Keplerian motion of the disk. 

Both NBO and HBO are observed on numerous occa- 
sions, but much importance has been given to the timing 
domain analysis (Cir X-1, Soleri et al. 2009; GX 340-^0, 
Penninx et al. 1991; GX 17-^2, Wijnands et al. 1996; 
Homan et al. 2002; Cyg X-2, Wijnands et al. 1997; GX 
5-1, Dotani 1988). Yu (2007) found HBO and NBO along 
with twin kHz QPOs in the source Sco X-1 and suggested 
that NBO resulted from the outward disk movement but 
there was no systematic spectral study to support this 
idea. In this Letter, we present the timing and spectral 
analysis of a RXTE observation of GX 5-1 during which a 
smooth transition of HBO-NBO-HBO is observed. The 
spectral study indicates that the inner disk edge moves 
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out radially or it gets thickened or puffed-up and the 
properties of the disk are substantially modified when 
compared to the properties of the disk with HBO. 

GX 5-1 is the second brightest Z source (Bradt et al. 
1968), located at a distance of 9.0±2.7 kpc (Christian & 
Swank 1997) with a luminosity in the range of 6.0 - 7.6 
X 10^*^ erg s-i (1 - 30 keV) (Jackson et al. 2009), and it 
resembles Cyg X-2 due to its HB nature (Kuulkers et al. 
1994). This source exhibits HBO and NBO along with 
kHz QPOs (van der Klis et al. 1985; Lewin et al. 1992; 
Wijnands et al. 1998) but FBO has not been detected 
yet. The kHz QPOs arc not seen at NB and FB which 
are often observed in Sco X-1. The detection of radio 
and infrared emission provides a clue for the existence of 
a jet (Fender & Hendry 2000; Jonkcr et al. 2000). 

2. DATA ANALYSIS AND RESULTS 

RXTE observed the source GX 5-1 on 1997 July 28 

spanning two satellite orbits (Obsld. 20055-01-04-00). 
We used the generic binned data (i?_250ws_4M_0 — 35_(3) 
for obtaining the power density spectrum (PDS) and the 
16 s binned data (standard 2 data) from the Proportional 
Counter Unit 2 (PCU2) for spectral studies; PCU2 is the 
best calibrated among the PCUs (Jahoda et al. 2006). 
We added 0.5% systematic errors to the spectral data 
to account for the calibration uncertainties. During this 
observation GX 5-1 was in the upper NB and the in- 
tensity was nearly constant in the first orbit (14913 ± 
50 c/s) with a HBO at 49±1 Hz but the intensity var- 
ied considerably in the second orbit, showing a dip-like 
feature. 

The left panel in Figure 1 shows the light curve of the 
second orbit in the energy range 2-13 keV, which is 
divided into 8 segments (a to h) for a detailed analy- 
sis. The intensity decreases to a lower state (~ 14000 
c/s) in about 1000 s and subsequently recovers to the 
initial state (^ 15000 c/s), similar to the variation seen 
in the X-ray dips of LMXBs. The right panels in Fig- 
ure 1 show the PDS in different energy bands, where 
the binned data are available in four energy bands i.e. 
1.94 - 3.70 keV, 4.05 - 5.82 keV, 6.18 8.68 kcV, 9.03 
12.99 kcV and the event mode for the 13.36 118.00 keV 
band. A NBO (6.0±0.7 Hz with a rms 2.00±0.22%) and 
a HBO (51. Oil. with a rms 2.28±0.30%) were found in 
the energy band 1.94 - 12.99 keV. The 6 Hz NBO and 51 
Hz HBO were not present in the lower (< 5.82 keV) and 
higher (> 13.36 keV) energy bands and were prominent in 
the intermediate energy bands i.e. 6.18 - 8.68 keV (7.1(t 
for HBO and 7.5cr for NBO) and 9.03 - 12.99 keV (3.2cr 
for HBO and 6.5cr for NBO). The HBO rms increased 
from 2.10±0.30% (6.18 - 8.68 keV) to 4.94±0.76% (9.03 
- 12.99 keV) and NBO rms increased from 1.92±0.23% 
to 4.07±0.49%. 

Figure 2 displays the best-fit PDS for the phases of 
ingress (or phase I: a+b+c), dip (II: d+e+f), and egress 
(III: g-l-h). In phase I the HBO was detected at a signif- 
icance level of 6.4cr (rms 2.58±0.31%) and its strength 
decreased in phase II (rms 1.55±0.35%) along with an 
appearance of NBO at 6 Hz (4.5a-; rms 1.92±0.42) and 
later HBO was observed (4.2(t) in phase HI. It is evident 
that the 6 Hz NBO is not present in phases I & HI. To 
study the evolution of PDS during the intensity variation, 
we obtained the PDS of each segment which is shown in 
Figure 3, along with the rms value and error. It is clear 



that the disappearance of HBO is almost balanced by the 
simultaneous appearance of NBO indicating a common 
physical origin behind the observed phenomenon. 

Both the eastern (a multi-color disk black body model 
plus a Comptonization model - diskbb-l-CompTT) and 
the western (a simple black body model, bbody, plus 
CompTT) approaches were used to fit the spectra of the 
three phases after excluding various other possible mod- 
els such as a single power-law, CompTT, diskbb-|-power- 
law, bbody-|-power-law, etc. The hydrogen equivalent 
column density was fixed at 6.6 x 10^^ cm~^ (Jackson 
et al. 2009). We confirmed that the allowance of the 
column density to be free does not improve the fit signif- 
icantly, indicating that the observed intensity variation is 
not related to absorption by dense blobs intervening the 
line of sight, as often observed in dipping X-ray sources. 
The top-left panels of Figure 4 show the best-fit spectra 
by the diskbb-|-CompTT model and the top-right pan- 
els represent the residuals when we apply the best-fit 
parameters of the phase I spectrum to the spectra from 
other phases. The residuals in the phase II spectrum sug- 
gest that the accretion disk properties have significantly 
changed at this phase. The best-fit spectral parameters 
of the three phases are given in Table 1 for both the 
adopted approaches. The source luminosity (in 2-10 
keV) varied from 4.54 x 10^^ erg s~^ (phase I) to 4.44 x 
10^® erg s"^ (phase II). 

It is difficult to meaningfully quantify the variations 
of spectral parameters because of the poor data statis- 
tics of the phase spectra, but residuals clearly indicate 
that there is a spectral diflference between phase I and 
II. To discern the difference, we investigated the mini- 
mum set of parameters which are sufficient and signifi- 
cant to justify the spectral change between the phases. 
We first fitted the phase I and II spectra (using the east- 
ern model) simultaneously with the phase II parameters 
tied to the best-fit parameters obtained from the phase I 
spectrum. The high value of x^/dof (==2998/104) clearly 
indicated that the spectrum has changed during the in- 
tensity variation (Table 1). We then found a decrease 
in the value of x^/dof as diskbb normalization (Ndisk) 
and CompTT normalization were allowed to be free, but 
the values of reduced were still unacceptable, indicat- 
ing the requirement of additional parameter variation. A 
substantial improvement of the fit was achieved by allow- 
ing the inner disk temperature (kTj„) to be free. After 
this step, any more significant improvement was not ob- 
tained even though all the other parameters were allowed 
to be free. The F-test values and probabilities indicate 
that the minimum set of spectral parameters needed to 
vary arc the normalizations of diskbb and CompTT and 
the inner disk temperature (Table 1). The simultane- 
ous fit results for all the three phases are summarized 
in Table 1. The analysis of contour plot shows that the 
disk normalization which is a close measure of inner disk 
radius {^disk oc R?„) varied significantly (at >99% con- 
fidence level, bottom- left panel of Figure 4). We applied 
the same method to the phase II and III spectra and 
found similar results (Table 1 and bottom-right panel of 
Figure 4). The study strongly suggests that the inner 
disk radius varied during the simultaneous evolution of 
HBO and NBO. 

3. DISCUSSION AND CONCLUSION 
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Wc have presented a detailed systematic study of the 
RXTE observation of GX 5-1 wlien the power of NBO 
and HBO were varying simultaneously (Figure 2 and Fig- 
ure 3). Yu (2007) found NBO, HBO, and twin kHz QPOs 
simultaneously in the source Sco X-1 and suggested that 
the generation of NBO is due to the radial outward move- 
ment of disk but this conclusion was not siipported by 
any spectral signatures. In this Lelter, we confirm such 
a disk movement during the evolution of HBO-NBO- 
HBO transition and propose that the associated inner 
disk edge has either significantly moved or puffed-up. 

We found that the ratio of NBO to HBO rms is close 
to 2 (panel d in Figure 3) similar to the value of Sco 
X-1 (Yu 2007) and this ratio drops to 0.4 in panel g 
(Figure 3). The absence of HBO and NBO in low en- 
ergy bands, relatively high amplitude of NBO in 9.03 
- 12.99 keV (Figure 1) and smooth transition of HBO- 
NBO-HBO (Figure 3) indicate that the QPOs are aris- 
ing from a hot region in the inner accretion disk. The 
simultaneous spectral fit result suggests that the inner 
disk edge moves out when NBO is dominant and other- 
wise it moves in (Figure 4 and Table 1). The observed 
spectral variations can be explained in the framework of 
the radiation-hydrodynamical model: the disk material 
leaves the inner edge and falls on to the neutron star sur- 
face covered with a hot corona. The outward radiative 
pressure halts or pushes back the inner disk edge, which 
can cause the relative increase of the observed disk nor- 
malization. Alpar et al. (1992) pointed out that the 
rotation of a thick disk in the inner region could cause 
the observed NBOs. In this scenario, however, the twin 
kHz QPOs should not be seen because of the slow dynam- 
ical movement but are often observed, e.g., in Sco X-1 
(where kHz QPO is also observed at normal and flaring 
branches). From the above discussion, we suggest that 
the oscillation frequency at the inner edge (at 51 Hz) is 
abated by the outward radiation pressure and results in a 
6 Hz NBO. Since the vertical height and surface density 
of the disk are expected to increase during the process, 
the observed increase in the disk normalization could be 
explained. 

In the two oscillator model, the perpendicular mode 
(vl) of oscillating hot blob in magnetosphere generates 
the HBO (Titarchuk et al. 1999) and if the mass of the 
blob increases, frequency of that mode would shift to a 
lower frequency not much afii'ecting the X-ray continuum. 
We speculate that whatever may be the oscillating mech- 
anism the observed smooth transition of HBO to NBO 
could be due to an increase in the participating mass in 
the oscillations. Our spectral results suggest that when 
the disk inner front radius which is the probable site for 
the transition boundary is increased, it gets thickened 
or puffed-up due both to the OTitward radiation pressure 
and to the accumulation of matter coming from the outer 
region of the disk. In this scenario, it is quite possible 
that the structure of the blob would be altered signif- 
icantly and the quasi-coherent perpendicular motion of 
the blob might have disturbed/randomized to produce 
the broad 6 Hz NBO instead of 51 Hz HBO. The other 
possible scenario could be that the Keplerian motion of 
the disk became super-Keplerian which affects the os- 
cillations (Titarchuk et al. 1998) in the; vertical struc- 
ture of the disk and have respective dependency on the 
oscillation frequencies {vy (viscous oscillations) oc 



(thickness) and z/t (break frequency) cx L~^, Titarchuk 
& Osherovich 1999). The NBOs arc always found to be 
< 20 Hz which implies that these flows would steeply 
affect these oscillations. 

Casella et al. (2005) associated HBO, NBO and FBO 
to C, B and A type QPOs in black hole binaries (BHBs). 
In general type B, A QPOs and their transitions are of- 
ten observed during the very high (VH) / steep power- 
law (SPL) state (McClintock & Remillard 2004) or in- 
termediate state (IM) (Belloni 2010). In this state, the 
accretion disk is considered to contain a compact corona 
(a low electron temperature with a high optical depth) 
along with a marginally truncated high temperature Ke- 
plerian disk (Done & Kubota 2006; Sriram et al. 2007, 
2010). A broad QPO at 6 Hz is often observed in a few 
BHBs e.g. XTE J1550-564 (Homan et al. 2001), XTE 
J1859-F226 (Casella et al. 2004), GX 339-4 (NespoH et 
al. 2003; Belloni et al. 2005). Similar transient variation 
from type B to A QPO was also observed in H1743-322 
(Homan et al. 2005) and GRS 1915-^105 (Soleri et al. 
2008). It was found that transient changes in QPOs are 
often associated with source intensity variations. The 6 
Hz NBO is observed only during a characteristic transi- 
tion in both Z-type sources and BHBs, which indicates 
that a common physical process is occurring at the inner 
region of the disk in both the sources. 

In a few atoll sources, similar QPOs varying between 5 
and 14 Hz with a relatively high coherence were reported 
(Wijnands et al. 1999; Wijnands & van der KUs 1999; 
Belloni et al. 2004). These detections and a low mass 
accretion rate in atoll sources pose a distinct challenge 
to the models explaining the 6 Hz NBO in Z sources. 
It was found that in atoll sources and BHBs these fea- 
tures are produced at the peak of their outburst where 
the inner disk edge is most probably located close to the 
compact object. We propose that the mass accretion 
rate is not the important criterion to produce NBO fea- 
tures, instead it is the location of the inner disk which 
drives the NBO. If the disk edge is close to the compact 
object, the chances of getting affected by the outward 
radiation pressure is high which would generate or dis- 
turb the oscillation in the inner region of the disk. Sim- 
ilar studies concentrating on the correlated spectral and 
temporal variation in LMXBs would be helpful to get a 
clear picture of the accretion disk during the presence of 
NBO/NBO-like features. 

We thank the anonymous referee for the useful com- 
ments. This research has made use of data ob- 
tained through HEASARC online service, provided by 
NASA/GSFC, in support of the NASA High Energy As- 
trophysics Programs. 
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Fig. 1.— Left: The RXTE/PCA 2-13 keV light curve (10 s 
binned) of GX 5-1 for Obsld. 20055-01-04^00. For detailed timing 
and spectral analysis, the light curve is divided into 8 segments (a 
to h). Right: The power density spectrum (PDS) of the light curve 
in different energy bands. The PDS in 1.94 - 12.99 keV is fitted 
by a model of a power-law with two Lorentzians. 
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Fig. 2. — The best-fit (dashed hnes) PDS for the ingress (or pliase 
I: a+b+c), the dip (II: d+c+f), and the egress (III: g+h) phase of 
the light curve. The Poisson level is not subtracted. The HBO and 
NBO rms and the errors (in brackets) are shown. 
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Fig. 3. — The PDSs of the 8 segments showing the temporal 
evolution of the coupled appearanee and disappearance of HBO 
(H) and NBO (N). The segment label and rms powers in QPOs 
(errors given in brackets) are given in each corresponding panel. 
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Fig. 4. — Top (left); The unfolded spectra for phases I (top), 
II (middle), and III (bottom), usmg the diskbb+CompTT model 
are shown. Top (right): Residuals resulting from modeling the 
spectra of all phases using the best-fit parameters for phase I spec- 
trum. Bottom (left and right): Significance contour plots between 
diskbb normalization (!^diak) of phase I and phase II, phase III and 
phase II obtained during the simultaneous fit procedure (see text) . 
Contours show the 68%, 90%, 99% confidence levels. Plots clearly 
suggests that ^disk is different for the three phases at more than 
99% confidence level. 
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TABLE 1 

Best-fit spectral parameters for the three phase spectra. 
The c;orresponding subscripts dbb k. bb are diskbb & bbody 
model. Errors are quoted at a 90% confidence level. The 
unabsorbed flux in 3.0 - 25.0 keV is given in units of 10~* 
ergs cm~2 s^i. The simultaneous fit results between 
phases I and II and phases II and III are shown. The 
best-fit parameter values from the simultaneous fit to 

THE three phase SPECTRA ARE ALSO LISTED. 



Jr'araiiit'lt'rs 


1 




11 






ill 








Ildbb 


Il66 


Illdtt 


III66 


kTir, (keV)'' 


2.13±0.05 


- 


2.00±0.04 


- 


2.17±0.05 






45±5 


- 


54±6 


- 


47±5 




knt {keVy 


- 


1.51±0.06 


- 


1.43±0.07 


- 


1.55±0.03 






0.079±0.016 




0.08±0.015 




0.10±0.002 


kTe'^'(kcV) 


3.21±0.05 


2.92±0.15 


3.13±0.05 


2.86±0.14 


3.53±0.07 


3.14±0.01 


T-' 


8.17±1.16 


9.0d±0.68 


8.21±1.24 


8.98±1.30 


7.33±0.85 


8.31±0.80 


Total FluX(di^fcf,fa+CompTT) 


3.63 




3.53 




3.63 




Total FluX(^bbody+CompTT) 




3.74 




3.57 




3.75 


X /dot 


44/ 46 


39/46 


43/ 46 


44/ 46 


36/46 


37/46 




Simultaneous Fit Result 














Phase I and II spectra 












Parameters 


5 

X 


dof 


F-fatat Value 


Probability 






All tied 


2998 


104 










Normalization's {Ndisk, ^CompTT) free 


899 


100 


58 


2.52x10 






kTin along with Norms free 


105 


98 


370 


2.01x10-*'* 






kTg, kTj^ along with Norms free 


92 


96 


6.7 


1.72x10-^ 






All free 


91 


92 


0.25 


0.91 








Simultaneous Fit Result 














Phase II and III spectra 












All tied 


2993 


104 










Normalization's (N^isk, ^CompTT) free 


778 


100 


71 


2.21x10"^** 






kTi„ along with Norms free 


100 


98 


332 


2.19x10"'"' 






kTe, kT^Ti along with Norms free 


86 


96 


7.8 


7.17x10"^ 






All free 


85 


92 


0.27 


0.89 








Simultaneous Fit Result 














Phase I, II, and III spectra 












Parameters 


I 


II 


III 








kTi^ 


2.13± 0.01 


1.96±0.01 


2.12±0.01 










44± 2 


66±3 


45±2 








CompTT Normalization 


8.95±3.28 


6.13±2.86 


7.54±3.29 









■^Inncr disk temperature of the diskbb model. 
''Normalization of the diskbb model. 

'^Temperature of the single temperature black body model. 

"^Normalization of the black body model. 

"^Electron temperature. 

^Optical depth of the Compton cloud. 
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